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1.0  Optoacoustic  Studies  of  Liquids 

1.1  Introduction 

Investigations  of  optoacoustic  generation  in  liquids  are  well  documented.  1  >2,3  Most  of 
these  studies  were  carried  out  using  low  optical  energies  and  small  optical  absorption 
coefficients.  A  small  optical  absorption  coefficient  makes  the  symmetry  of  the  problem  simple; 
the  mathematics  tractable.  In  the  low  optical  energy  regime,  the  acoustic  generation  mechanism 
is  thermal  expansion.  This  is  the  basis  for  linear  optoacoustic  theory.  Nonlinear  effects  occur 
when  the  liquid  experiences  multiphoton  absorption  or  a  phase  transition. 

On  a  molecular  level,  the  optoacoustic  effect  can  be  broken  down  into  several  different 
steps.  The  process  of  converting  optical  energy  to  acoustical  energy  begins  with  the  absorption 
of  photons  by  the  molecules  of  the  fluid.  This  absorption  leaves  the  molecules  in  an  excited 
state.  In  the  present  case,  a  nitrogen  laser,  used  as  the  energy  source,  produces  electronic 
excitation  of  the  molecules.  Dissipation  of  this  energy  can  occur  through  several  mechanisms 
including;  photochemical  processes,  radiative  decay,  and  nonradiative  decay.  In  nonradiative 
decay,  the  electronic  excitation  energy  is  transferred  to  translational,  rotational  or  vibrational 
states  of  the  molecule.  Through  interactions  with  other  molecules  this  rotational  and  vibrational 
energy  can  be  converted  into  translational  energy  creating  a  temperature  rise  in  the  fluid.  In  our 
experiment,  following  each  pulse  from  the  laser,  the  fluid  is  allowed  to  cool  to  an  ambient 
temperature.  This  alternate  heating  and  cooling  creates  a  pressure  variation  which  is  observed 
as  acoustic  pulses. 

Present  optoacoustic  theories4>5,6  assume  that  the  energy  transfer  processes  discussed 
above  all  take  place  over  a  period  of  time  which  is  very  short  compared  to  other  factors  that 
effect  the  observed  pressure  profile;  such  as,  the  length  of  the  laser  pulse  and  its  diameter  as  it 
enters  the  fluid. 

The  process  of  energy  transfer  and  the  rate  at  which  it  occurs  varies  from  one  liquid  to 
another  because  the  energy  levels  and  rate  of  energy  transfer  depend  on  the  structure  of  the 
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molecules.  If  the  rate  of  production  of  translational  energy  is  slow  due  to  internal  relaxation 
processes,  the  shape  of  the  acoustic  pulse  should  be  governed  by  these  relaxation  processes. 
The  experimental  system  must  be  capable  of  resolving  very  short  acoustic  pulses  in  order  to 
resolve  events  occuring  on  this  time  scale.  Otherwise,  the  pressure  profile  is  dominated  by 
geometrical  factors  of  the  experiment. 

Carbon  disulfide  (CS2)  has  been  chosen  as  one  subject  of  study  because  previous 
ultrasonic  absorption  measurements7  have  already  shown  that  the  transfer  of  energy  between 
translational  and  vibrational  states  is  relatively  slow  (compared,  for  example,  to  H2O).  This 
transfer  of  energy  into  translation  can  be  described  by  a  relaxation  equation  which  gives  rise  to 
an  exponetial  decay  in  the  number  of  excited  molecules.  Since  some  of  the  energy  will  go 
directly  from  electronic  to  translation,  it  is  expected  that  the  optoacoustic  signal  in  CS2  will  be  a 
simple  superposition  of  such  an  exponential  decay  and  a  signal  such  as  occurs  in  liquids  where 
the  energy  transfer  takes  place  much  more  quickly.  By  extracting  this  exponential  decay 
process  from  the  raw  data,  we  can  determine  the  relaxation  time  that  characterizes  the  energy 
transfer  process. 

One  can  view  a  liquid  as  being  composed  of  dynamic  molecular  groups  referred  to  as 
clusters.7  The  average  number  of  molecules  per  cluster  differs  for  various  fluids.  Photons 
absorbed  by  these  clusters  cause  a  redistribution  or  reorientation  of  the  molecules  which  results 
in  a  change  in  density.  These  molecular  movements  result  in  an  outwardly  spreading  acoustic 
wave.  By  monitoring  characteristics  such  as  amplitude,  relaxation  time,  and  the  shape  of  the 
wave,  one  might  be  able  to  gain  insight  to  the  molecular  dynamics  of  the  system.  For  example, 
a  state  of  increasing  instability  is  expected  with  increasing  photon  density.  As  Np  (the  number 

of  photons  absorbed  per  unit  volume  in  one  relaxation  time)  surpasses  the  number  of  clusters 
(Nc),  clusters  begin  to  break  up  to  a  greater  extent  until  the  liquid  vaporizes.  As  the  system 
nears  this  regime,  quantities  characteristic  of  the  liquid  state  approach  those  of  the  gaseous 
state.  It  is  near  this  phase  transition  that  the  established  linear  theory  is  no  longer  applicable. 
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Experimental  Configuration 


A  block  description  of  the  experimental  technique  is  shown  in  Figure  1.1  This  technique 
utilizes  two  lasers.  A  1.25  MW  nitrogen  laser  emitting  800  ps  pulses  in  the  UV  (337  nm)  is 
used  as  the  optoacoustic  source.  A  3mW  continuous  He-Ne  laser  runs  parallel  to  the  excitation 
zone.  It  is  then  focused  on  a  photodiode  which  is  coupled  to  a  digitizing  oscilloscope.  The 
probe  laser  responds  to  refractive  index  gradients  which  are  produced  as  a  result  of  the  acoustic 
wave.  These  gradients  cause  small  deflections  of  the  probe  beam  which  are  detected  by  the 
photodiode.  The  observed  signal  is  proportional  to  the  pressure  gradient.  Comparison  to 
theory  can  be  achieved  by  either  integrating  the  observed  signal  or  by  differentiating  the 
predicted  pressure.  Both  techniques  have  been  used;  however,  a  comparison  to  the  absolute 
magnitude  of  the  acoustic  pressure  can  be  achieved  as  well  if  the  theory  is  compared  to  the 
output  of  a  PVDF  hydrophone  as  shown  in  Figure  1.2.  Given  the  inaccuracies  in  beam  size 
estimates,  the  agreement  between  theory  and  experiment  is  considered  acceptable. 

Energy  delivered  to  the  fluid  is  varied  by  placing  glass  slides  in  the  beam  path  of  the  UV 
laser.  N  number  of  glass  slides  (~  .23  mm  thick)  are  supported  by  two  small  pieces  of 
plexiglass.  The  centers  of  the  plexiglass  have  been  removed.  To  provide  for  accurate  energy 
measurements,  twenty  sets  of  attenuators  corresponding  to  a  different  number  of  glass  slides 
were  made.  The  energy  transmission  for  each  attenuation  was  measured  with  an  energy  meter. 
The  precision  of  the  energy  meter  is  +  lpJ.  Figure  1.3  shows  measured  energy  delivered  to 
the  test  cell  as  a  function  of  the  number  of  glass  slides. 

The  energy  density  in  the  fluid  can  also  be  increased  by  increasing  the  optical  absorption 
coefficient  a.  This  is  achieved  by  either  adding  dye  to  the  test  liquid  or  by  choosing  a  liquid 
that  strongly  absorbs  UV  radiation. 

1.2  Initial  Studies  Near  the  Non-Linear  Regime 
In  optoacoustical  studies  on  the  molecular  level,  among  the  parameters  that  must  be 
determined  are  the  number  of  fluid  molecules  being  excited  and  the  number  of  photons  being 
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absorbed  per  unit  time.  The  number  of  photons  is  easily  obtained  by  measuring  the  total 
energy  of  the  pulse  with  an  energy  meter  and  dividing  by  the  energy  per  photon.  The  number 
of  molecules  is  obtained  by  first  measuring  the  optical  absorption  coefficient  a  of  the  liquid. 
The  optical  penetration  depth  is  then  the  reciprocal  of  a.  Knowing  the  laser  beam  dimensions 
and  assuming  exponential  decay,  one  can  then  find  the  excitation  volume. 

The  optical  absorption  coefficient  of  pure  CS2  at  X  =  337  nm  is  large.  To  measure  a, 
small  amounts  of  CSq  were  added  to  propanol  which  has  a  much  smaller  absorption 
coefficient.  The  energy  passing  entirely  through  the  test  cell,  measured  by  an  energy  meter,  is 
then 

E  =  E0e'[(l'f)a,  +  fa^x'a  where  .. 


al  =  absorption  coefficient  of  pure  propanol  and  the  test  cell 
otq  =  absorption  coefficient  of  pure  CS2 
a  =  absorption  of  the  test  cell 
f  =  fraction  of  CS2  in  propanol,  and 
x  =  length  of  test  cell  (4  cm). 

In  all  cases  of  interest,  f  and  a  are  very  small  so  oc^l-fjx  +  a  =  <xix.  The  quantity  ccix  can  be 

measured  by  first  measuring  the  energy,  Emax,  in  the  absence  of  the  test  cell  and  then  the 
energy,  E0,  when  the  test  cell  full  of  propanol  is  placed  in  the  beam  path.  We  then  have 

p  p  -0C  jX 

”  ^maxe  #  (12) 

For  our  system,  Emax  was  found  to  be  0.363  mJ  and  E0  was  found  to  be  0.340  mJ  so  alX  = 
0.0655.  We  now  can  write  Eq.  (1.1)  as 


t,  ^  -0.0655 

E  =  E0e  e 


-fa  2 


(1.3) 


The  total  volume  V  of  propanol  for  this  experiment  was  12  mL.  Drops  of  CSq  with  a 
volume  of  0.01  ml  were  added  to  the  propanol  so  that 
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f  N-(.Ol) 
12 


(1.4) 

where  N  is  the  number  of  drops.  Energy  transmission  was  measured  after  each  drop.  From 
this  data  a2  was  deduced  to  be  370  cm'1  which  implies  a  penetration  depth  for  pure  CS2  of 
27pm. 

The  excitation  volume  is  then  given  by 

V  =  xy  f  e  a^dz  where 

J°  (1.5) 

x  and  y  are  the  spatial  dimensions  of  the  rectangular  laser  pulse.  From  the  measured  values  of 

x,  y,  and  a2  we  obtain  B  =  2.22  x  10-6cm3. 

The  number  of  molecules  n,  in  the  excitation  volume  is 


n  = 


pVN 

W 


where 


p  =  density  of  CS2  =  1.263  g/cm^ 

N  =  Avagadro's  number 
W  =  molecular  weight  of  CS2 

From  these  values,  the  excitation  volume  contains  2.22  x  1016  CS2  molecules. 


(1.6) 


A  nonlinear  transition  has  been  observed  in  CS2  when  sufficient  energy  density  is 
applied.  We  now  refer  to  Figure  1.4  which  is  a  plot  of  acoustic  amplitude  vs  laser  pulse 
energy  in  CS2.  There  are  three  regimes  of  interest  on  this  curve  which  are  separated  by  two 
regions  of  constant  acoustic  amplitude.  The  first  extends  from  zero  to  75  pj.  The  boiling 
point  of  CS2  is  46°  C.  The  equilibrium  temperature  of  the  liquid  was  22°  C.  In  order  to 
achieve  the  boiling  point,  a  temperature  rise  of  24°  C  is  required.  Knowing  the  specific  heat 
and  the  mass  of  the  fluid,  one  can  determine  the  laser  pulse  energy  required  to  achieve  the 
boiling  point.  The  boiling  point,  bp,  is  indicated  near  this  first  plateau.  The  second  regime 
extends  from  90  |lJ  to  160  pJ.  The  slope  of  the  second  regime  increases  relative  to  the  first 
indicating  a  higher  optoacoustic  conversion  efficiency  probably  due  to  vaporization  of  the 
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liquid.  The  point  where  vaporization  bubbles  appeared  to  the  unaided  eye  is  indicated  by  vb. 
The  third  regime  extends  from  190  jlJ  to  Emax.  The  slope  changes  again  indicating  optical 
saturation  of  the  liquid. 

1.3  Opto-Acoustic  Observations  in  a  Relaxing  Fluid 

Ultimately  we  wish  to  use  the  optoacoustic  effect  to  investigate  relaxation  processes  in 
fluids.  For  most  fluids,  internal  reorientations  take  place  very  rapidly,  much  faster  than  the 
present  experimental  apparatus  can  resolve.  As  an  interim  step,  fluids  with  relatively  slow 
internal  relaxations  will  be  explored.  The  prototype  for  such  behavior  is  CS2.  Ultrasonic 
attenuation  measurements  have  shown  that  liquid  CS2  has  a  very  slow  internal  relaxation.8 

For  relaxation  measurements,  nonlinearity  is  undesirable  so  a  sufficient  number  of  glass 
slides  were  placed  between  the  excitation  laser  and  the  CS2  test  cell  to  insure  operation  in  the 
linear  regime.  Even  in  the  linear  regime,  the  large  absorption  coefficient  presents  difficulties. 
Specifically,  since  the  excitation  energy  is  all  absorbed  within  a  depth  not  much  different  in  size 
from  the  laser  beam  waist,  the  source  for  the  optoacoustic  signal  appears  more  spherical  than 
cylindrical.  The  more  complex  geometry  results  in  theoretical  expressions  which  are  much 
more  complex  than  those  used  by  previous  investigators.^ 

Experimental  Results 

Prior  to  launching  into  a  tedious  mathematical  derivation,  it  is  of  interest  to  examine  the 
general  shape  of  the  experimental  curves.  Figure  1.5  shows  the  magnitude  of  the  probe  beam 

O 

deflection  versus  time  for  propanol.  Hutchins  and  Tam  found  this  same  general  shape  for 
other  fluids  as  well.  In  this  case,  energy  is  transferred  to  the  fluid  more  rapidly  than  the 
experimental  apparatus  can  respond  so  the  waveform  is  dictated  by  laser  beam  geometries. 
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Also  shown  in  Figure  1.5  is  a  wave  form  observed  under  similar  conditions  for  CS2. 
Note  that  although  the  rise  time  is  comparable  to  that  for  propanol,  the  decay  is  much  slower 
and  the  waveform  spends  much  more  time  positive.  It  is  not  yet  clear  whether  this  behavior  is 
a  result  of  a  slowly  relaxing  internal  mode  or  the  much  larger  optical  absorption  in  CS2  which 
leads  to  a  more  spherical  source  term. 


Relaxation  Equations 

Assume  that  the  incoming  laser  beam  is  radially  Gaussian  with  a  beam  waist  a,  with  a 
Gaussian  time  dependence  characterized  by  tp,  and  is  absorbed  in  accordance  with  Beer's 
law, ^  then  the  intensity  I  can  be  written  as 


I(r,z,t) 


I„Ve 


t  exp 


2  2 

-r  t 

- az  — 

2  2 

a  t 


VI 


(1.7) 


The  energy  absorbed  per  unit  time  per  unit  volume  by  the  fluid  is 

w(r,z,t)  =  al(r,z,t)  ^  g) 

but  this  is  not  the  energy  which  appears  in  translation  per  unit  time.  This  fact  makes  our 
treatment  different  from  previous  ones.  We  will  assume  that  only  when  the  energy  appears  in 
translation  is  there  a  change  in  density.  In  fact,  an  increase  in  electronic  or  vibrational  states 
could  lead  to  density  changes  through  variations  in  bond  lengths  or  intermolecular  bonding,  but 
these  will  be  ignored  for  now. 
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For  simplicity,  consider  a  three  level  system  as  shown  below. 
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Following  the  notation  of  Bauer,  ^  we  introduce  a  set  of  progress  variables, 

•  * 

^a=n(kaxixj-kaxkx1) 

where  a  refers  to  an  energy  transfer  reaction,  i  - 1  refer  to  initial  and  final 
pairs,  and  ka  and  ka'  are  forward  and  reverse  reactions  respectively.  In 
change  in  fluid  translational  energy  is  given  by 

Ea=^aAEa  (1.10) 

where  AEa  is  the  energy  deficit  in  reaction  a. 

It  should  be  noted  at  this  point  that  in  a  fluid,  ternary  collisions  as  well  as  unimolecular 
reactions  might  be  important.  These  can  be  handled  with  pseudoreactions  of  the  above  type 
provided  we  are  always  operating  at  the  same  liquid  density. 

At  equilibrium, 

°  °  t/°°  °-n 

KaxixJ  "  Kctxkxl  ~  u  •  (-I  in 


t,a  such  that 

(1.9) 

states  of  collision 
this  notation,  the 


If  we  expand  near  equilibrium, 


y  1  ,  0  0  0  dkaxixj  —SlCojXjXj 

Wn  =  kaxixj  +—  TT~  dT+2, - dn 


3T 


dn 


m 
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O'  o  o  9kaxkx/  _  K,dkaXkXl 


o  o  O  0 r'K'M  ™  V 

-kaxkx/ - dT  ~2j 


dT 


dn. 


dnm. 


(1.12) 


Note  that  ka  depends  upon  T  not  nm  while  x's  depend  upon  nm, 

.  ,  '  -AE(/RT 

"koc  =  kae 


(1.13) 


and 


ni 

xi  =  — . 

1  n 


(1.14) 


We  can  then  write 


dkaxixj  jrr  ,  o  o^a 
- dT  =  kaXjXj - dT  ,  etc. 

dT  3T 


(1.15) 


so 


■  \ 

3lnk  dink 

uuix>.a  a 


dT  dT  / 


0  O  oY 
_X:X:  Is 


+  kCCXiX 


dlnXjXj  3lnxkx; 


“  \  dnm  dnm  y 


dn. 


(1.16) 


or 


Sjn^xf 


dT 


-Mkc/kaXT-To) 


+  2< 

m 


dlnxj  3lnXj  3lnxk  3lnx^ 


y  dnm  dnm  dnm  dnm 


(nm "  nm) 


(1.17) 


or 


ij »  -  on!  41-  w) (t  -  T”)+ A . .  1 . 


dT 


m  dnm  dnm 


(1.18) 


14 


Any  dependence  of  on  nm  is  explicit  except  for  m  =  i.  So 

0 


^  3lnnj  /  o  \ 

m  onm 


With  this  observation  included, 

AE 


-  k^xj  ^  (T .  t")  +  +  2LJSL .  5i^l  -  ^ 


0  0  0 
k 


RT 


ni  nj 


n,  nk 


We  are  interested  in  reactions  of  the  form 

k, 


a 


nij  +  mj  m/+  mj  +  AEa , 

that  is,  j  =  k,  so  the  process  takes  place  at  the  same  rate  no  matter  what  the  collision 
Since  our  reactions  may  be  unimolecular,  this  form  is  especially  desirable  since  now, 


,  af 

o  o  o  / 


RT 


A  o  o 

0  Ti;  ”  n :  n/-n/ 

T-T  +  1  1  11 


We  are  dealing  with  three  reactions, 


and 


o  o  j  i 

RT2 


J.  .  O  oj  AE2 

S2=k20*2x< 

RT" 


(t-tJ 


+ 


0  oj  AE12 
S3  =  k21x2xl 

RT" 


(t  -  T°)  +  — 


ni 

n,  j 

0 

nl  ■  nl 

r 

O  O 
£ 

O 

a 

nl 

n0 

0 

n2  '  n2 

o3 

1 

3 

0  0 

n2 

no 

0 

n2  '  n2 

nl  ■  nl 

n2 


(1-19) 


(1.20) 


partners. 


(1.21) 


(1.22a) 


(1.22b) 


(1.22c) 
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From  these  equations,  one  can  solve  for  relaxation  times,  X ^  and  %2  and  relaxation  strengths, 
Ai  and  A2  which  completely  describe  the  transfer  of  energy  to  translation.  Actually,  one  other 
equation  is  needed, 

n  =  no  +  ni  +  n2  (1.23) 

and  Aj  and  A2  depend  upon  the  excitation  source.  Solutions  of  this  equation  will  be  left  for 
later. 


Wave  Equation 

The  inhomogenous  wave  equation  for  acoustic  pressure,  ignoring  electrostriction,  can  be 
written  [according  to  Eq.  (15)  of  Ref.  2]  as 


2 

2  _  2 
\c  dt 


\ 

2 

-  V 

/ 


p  = 


a(3  d 
Cp  dt 


I(r,z,t) 


(1.24) 


where  c  is  the  speed  of  sound,  a  is  the  optical  attenuation  coefficient,  [3  is  the  expansion 
coefficient,  Cp  is  the  specific  heat  and  I  is  the  source  intensity.  We  can  make  use  of  this 
equation  if  we  recognize  that  al  is  the  rate  of  energy  deposition  per  unit  volume  which,  in  our 
case  (ijaser  pU|se  «  ^relaxation)  controlled  by  energy  transfer  (relaxation). 


Since  tp  (same  as  Tjaser  pUise)  is  very  small,  we  will  rewrite  the  equation  for  intensity  as 


I(r,z,t)  = 


-t/T 

e 


+  A 


-t/T 


2/T7 


2  2 

-r  /a  -az 
e  e 


(1.25) 


This  is  no  longer  the  laser  intensity  but,  rather,  is  the  laser  intensity  modified  for  the  time 
required  for  energy  transfer.  Normalizing, 
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Ep(energy  in  pulse)  =  |l  (r',z',t')  r'dr'dO'dt 


(1.26) 


or 


f  °°  _ f’/T  _T-'  fn 

Ep=|  (A1/Xe  )  +  e  r'dr'dG'dt' +  like  terms  for x2. 

J  n 


(1.27) 


Integrating  over  all  time  and  space, 

Ep  =  7ta  (Aj  +  A2}  =$  A2  =  — 2- .  Aj. 

Tta 


(1.28) 


It  is  sometimes  convenient  to  introduce  a  scalar  potential,  (j>,  at  this  point  such  that 


and 


2 

1  a 

2  ..  2 

\c  3t 


<t>  =  I(r,z,t) 


then 
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ap  3$ 

cp  3t 


1  a2-v2 
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(N 
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c  dt 

L 

^1 


T1A2 


-t/T 

e 


4-  e 


-t/T  * 


2  2 

-r  /a  -az 

e  e 


(1.29) 


(1.30) 


(1.31) 


This  inhomogeneous  wave  equation  can  be  solved  using  Green's  functions.  The  general 
form  is 


2 

V  \j/- 


2 

3  \)/ 


3t 


-47tq  (r,t) 


(1.32) 


17 


in  which  case 


q(r0,  t-R/c) 


V(M)--7-fdVo:£1V 

471-'  ^ 


(1.33) 


where 


R  = 


(x-x0)2  +  (y-y0)2  +  (z-Zof 


1/2 


(1.34) 


If  we  transform  to  cylindrical  coordinates, 


x  =  r  cos  0  y  =  rsin0  z  =  z 


R  =  Vr2  +  r2  -  2rr0  cos(0-0o)  +  (z-z0)2  dV0  =  rodrod0 0dz0  so 


,  ,  47tx2 

<Kr,t)  =  —  dV0- - 

4k  J 


AlT2  -(t-R/cJ/X,  -(t-R/cJ/x ; 
C  “1*  c 


A2X! 


e^V“z° 


R 


or 


4>(r,t)  = 


T?(47t) 


f  "l 

J  dV0  — 


X2  -(t-R/c)/X]  -(t-R/c)/x; 

*  6  -I-  C 

t.  A2 


e"H"- 


R 


(1.35) 


and 
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2 

-AiX2  -(t-R/c)/x,  -(t-R/c)/x2 

-r%:  -cxz0 

cc(3  A2 

f  dV  - 

2 

A2  xl 

2  2  . 

(4k)  Cp 

1  u  v  0 

R 

We  measure  dn(r,t)/dr  or,  if  we  assume  n  a  p  a  p,  then  3n(r,t)/9r  is  proportional  to 


dp 

dr 


gp  a2 

2  2 
(47C)  Cp  %2 


I 


dV, 


2 

A^2 


A2tj 


-(t-R/c)^!  -(t-R/cJ/T! 

c  c 


cXjR 


-(t-R/c)/x2  -(t-R/c)/x2  ] 

c  c 


R~ 


R 


r  -  r0  cosCG-G^J  e’r°/ae  az°. 


Current  Status 

At  this  time,  we  have  observed  an  optoacoustic  signal  in  CS2.  The  optoacoustic  signal  is 
quite  different  from  that  observed  in  propanol.  Mathematical  equations  taking  into  account 
relaxation  and  strong  optical  absorption  have  been  derived  and  programmed.  We  are  now  in 
the  process  of  comparing  experimental  results  to  these  calculations.  This  investigation  will 
comprise  an  M.S.  thesis  by  Mr.  Charles  Thompson  with  an  expected  completion  date  of  15 
December  1987. 
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2.0  Ultrasonic  Attenuation  in  Polymer  /Water  Mixtures 


During  the  period  1  May,  1985  -  30  September,  1985  in  cooperation  with  Professor  John 
Jacobus  at  Tulane  and  with  support  by  the  Naval  Oceanographic  Research  and  Development 
Activity  ultrasonic  attenuation  was  measured  in  a  variety  of  aqueous  polymer  solutions.1  This 
initial  study  identified  three  polymers  as  having  significant  attenuation  in  excess  of  that  for  pure 
water.  These  were  Polyethylene  Oxide  (MW  =  14,000),  Polycoethylene  Glycol  (MW  = 
15,000),  and  Hexadecyl  trimethyl  ammonium  bromide.  All  showed  excess  attenuation  near 
2.5  dB/m  at  some  frequencies,  but  the  14K  PEO  showed  the  maximum  attenuation  at  630 
KHz. 

The  ultrasonic  apparatus  used  in  the  preliminary  study  was  restricted  to  a  range  0.63 
MHz  to  2.5  MHz  with  an  absolute  accuracy  of  approximately  1.0  dB/m.  Though  this  accuracy 
was  acceptable  at  the  time  since  our  goal  was  to  screen  solutions  for  those  having  an  excess 
attenuation  of  10  dB/m,  the  frequency  range  was  not  sufficient  to  identify  mechanisms.  These 
preliminary  measurements  illustrated  that  a  solution  with  good  attenuation  at  high  frequency 
(CTAB  at  1.6  MHZ)  might  show  almost  no  attenuation  at  lower  frequencies  (630  KHz)  which 
is  suggestive  of  a  relaxation  process.  The  reverse  is  also  true  (14K  PEO  for  example). 

In  the  past  few  months,  an  ultrasonic  apparatus  designed  to  operate  at  lower  frequencies 
was  assembled  and  tested.  This  apparatus  was  originally  proposed  for  the  initial  work  but 
delays  in  funding  prompted  us  to  go  with  an  available  instrument  instead  of  the  optimum 
system. 

2.1  Ultrasonic  Apparatus 

The  ultrasonic  apparatus  to  be  used  for  chemically-induced  ultrasound  attenuation 
measurements  (Figures  2. 1-2.3)  consists  of  sending  and  receiving  transducers  separated  by  a 
variable  distance,  Ax.  The  transmitting  transducer  is  driven  by  an  amplified  tone  burst 
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Figure  2.1 
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RECEIVING  TRANSDUCER  ASSEMBLY 


Figure  2.2 
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Figure  2.3 


generated  by  a  variable  frequency  oscillator  (tone  burst  generator)  -  amplifier  combination.  The 
burst  of  sine  waves  (approx.  5  cycles  at  93  KHz,  100  cycles  at  1  MHz)  can  be  detected  by  the 
receiving  transducer. 

The  received  signal  is  preamplified,  amplified  and  filtered  before  it  is  fed  into  an 
oscilloscope  where  the  amplitude  is  measured. 

2.2  Experimental  Results 

The  earlier  measurements  indicate  that  a  20  gm/liter  solution  of  14K  PEO  in  water  will 
provide  an  attenuation  of  2.8  ±  0.5  dB/m  at  630  KHz.  The  new  instrument,  then,  should  be 
capable  of  measuring  attenuation  to  an  accuracy  of  ±  0.5  dB/m  from  below  100  KHz  to  at  least 
630  KHz  (to  provide  overlap  with  previous  measurements).  In  order  to  test  accuracy,  a  series 
of  measurements  was  conducted  in  distilled  water.  The  results  are  given  in  Table  2.1.  As  can 
be  seen,  the  standard  deviation  (a)  is  an  order  of  magnitude  lower  than  the  maximum 
acceptable  value  demonstrating  that  this  system  is  capable  of  the  required  accuracy. 


Measured  Absorption  in  Distilled  Water 
(Average  of  Five  Measurements) 


Freq.  (KHz) 

a(dB/m) 

a(dB/m) 

91 

0.028 

0.038 

172 

0.012 

0.008 

508 

0.024 

0.052 

851 

0.014 

0.018 

1191 

0.004 

0.016 

Table  2.1 


2.3  Status  of  Research 

Measurements  are  now  underway  in  PEO.  All  measurements  and  analyses  are  scheduled 
to  be  competed  by  the  end  of  this  year,  when  this  sub-task  will  be  terminated. 

2.4  References 

1.  "Chemically  Induced  Acoustic  Attenuation,"  by  Brenda  Little,  Bert  Green,  O.  John 
Jacobus,  and  Henry  Bass,  Final  Report,  October,  1985. 
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3.0  The  Propagation  of  Sound  in  Vibrationally  Excited  Gases  (SACER) 

3.1  Summary  of  Past  Work 

Bauer  and  Bass  predicted  in  1973  the  amplification  of  sound  in  a  gas  with  an  over¬ 
population  of  vibrationally  excited  states  and  named  the  effect  "SACER"  (sound  amplification 
from  controlled  excitation  reactions)1.  Calculations  were  made  in  1984  2  of  the  magnitude  of 
the  effect  in  a  number  of  gases  where  the  relaxation  was  slow  enough  to  allow  the  propagation 
of  a  sound  wave  in  the  excited  gas.  During  the  past  year,  the  results  of  measurements 
demonstrating  the  phenomena  in  N2/H2  mixtures  were  published,3  and  initial  measurements  in 
N2/He  and  N2/CH4  mixtures  were  reported.^  Measurements  have  also  now  been  made  on 
N2/H2O  mixtures.  Because  the  relaxation  times  in  these  mixtures  vary  widely,  both  in 
magnitude  and  temperature  dependence,  this  series  of  measurements  serves  as  a  strenuous  test 
of  the  theory. 

Initial  measurements  have  also  now  been  made  on  CO/H2  mixtures.  For  these 
mixtures  an  effort  will  be  made  to  monitor  the  decay  of  the  vibrational  temperature  of  the  CO 
by  the  infrared  emitted  by  the  excited  gas. 

3.2  Experimental  Method 

The  experimental  equipment  is  diagrammed  in  ref.  3.  An  abrupt  electrical  discharge 
through  the  gas  contained  in  a  tube  2.54  cm  in  diameter  and  about  30  cm  long  leaves  the  gas  in 
a  vibrationally  excited  state  and  produces  a  resonant  acoustical  vibration  in  the  tube.  The  sound 
velocity,  and  hence  the  translational  temperature  in  the  gas,  is  determined  by  the  change  in 
frequency  of  the  oscillations  and  the  sound  absorption  (or  amplification)  is  determined  from  the 
change  in  amplitude  of  the  oscillations.  Samples  of  the  decaying  acoustical  oscillations  are 
shown  in  Fig.  3.1. 

In  all  of  the  mixtures  except  N2/H2O,  the  mixtures  were  prepared  and  stored  in  a 
supplemental  tank.  For  N2/H2O  mixtures  the  desired  concentration  was  obtained  by  bubbling 


26 


Filtered  Microphone  Response  (Arbitrary  Units) 


1 

G 

1 

V  40  Torr 

7.3  kJ/mol 

— L —  J - 1 _ 1 _ 1 _ 

o  60 


o 


60 


60 


TIME  (ms) 


Figure  3.1  Decay  of  acoustical  oscillation  in  vibrationally  excited  mixtures  of 
N2  with  H2>  CH^  and  H20.  The  concentrations,  gas  pressures  and  initial  excess 

energy  in  vibration  is  shown  on  each  figure. 
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the  N2  through  an  ice  water  trap  at  the  pressure  needed  for  the  specific  concentration.  The 
mixture  was  allowed  to  flow  through  the  sound  tube  until  water  vapor  concentration  had 
equilibrated  with  the  tube  wall. 

The  measured  sound  velocity  is  corrected  for  the  tube  wall  loss  and  then  used  to 
determine  the  variation  of  the  translational  temperature  with  time  following  the  electrical 
discharge.  Some  of  these  temperature  vs.  time  curves  are  shown  plotted  as  points  in  Fig.  3.2. 

In  order  to  analyze  the  results  it  is  desirable  to  know  the  initial  vibrational  temperature 
and  its  decay  rate.  In  N2  mixtures,  the  initial  Tv  is  obtained  by  measuring  energy  deposited  in 
the  gas  by  the  discharge  and  assuming  that  practically  all  of  this  energy  goes  initially  into 
molecular  vibration.  For  this  purpose  the  voltage  pulse  across  the  tube  and  the  current  pulse 
through  the  tube  are  measured  and  stored  by  a  digital  oscilloscope.  The  energy  is  determined 
by  integrating  the  product  of  these  curves.  The  discharge  lasts  approximately  250  ^seconds, 
and  the  storage  scope  samples  at  a  0.5  |isec  rate.  The  decay  of  the  vibrational  temperature  is 
inferred  from  the  measured  time  dependence  of  the  translational  temperature  using  the 
relaxation  model  discussed  below. 

For  the  case  of  CO  mixtures  it  should  be  possible  to  get  a  direct  measure  of  the 
vibrational  temperature  as  a  function  of  time  from  a  measure  of  the  infrared  radiated  by  the 
excited  gas.  For  this  purpose  a  CaF2  window  has  been  mounted  on  the  end  of  the  discharge 
tube  and  the  radiation  coming  out  of  the  tube  focused  on  a  photo  detector.  To  date  the  decay  of 
the  infrared  has  been  observed  but  the  technique  for  extracting  an  absolute  temperature 
measurement  from  it  has  not  been  perfected.  Difficulty  has  also  been  encountered  in  getting 
enough  energy  into  the  CO  gas  from  the  electrical  discharge.  The  CO  seems  to  have  greater 
conductivity  and  the  present  transformer  is  not  capable  of  delivering  enough  current.  Though 
theoretical  calculations  indicate  that  CO  should  be  similar  to  N2,  it  may  be  that  more  of  the 
discharge  energy  is  going  into  ionization  in  the  CO. 
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Figure  3.2  Variation  of  the  translational  temperature  in  the  gas  mixtures 
shown  in  Fig.  3.1.  The  points  were  calculated  from  the  sound  frequencies 
obtained  from  the  pulses  in  Fig.  3.1.  The  curves  were  calculated  as  explained 
in  the  text. 
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During  the  past  summer  a  four  channel  digital  scope  has  been  acquired  which  allows 
the  simultaneous  storage  of  the  discharge  voltage,  the  discharge  current,  the  infrared  decay, 
and  the  acoustic  decay. 

3.3  Analysis  of  the  Results 

The  experimental  translational  temperature  (the  points  in  Fig.  3.2)  have  been  compared 
to  curves  calculated  from  the  following  relaxation  model: 

X, 


vVC 


-v 

i 


-»  wall 


->  wall  . 


Ltc 

Here  the  vibrational  energy  is  assumed  to  decay  with  the  passage  of  energy  into  translation  and 
to  the  tube  wall  with  characteristic  decay  times  of  xy  and  xvc  respectively.  The  translational 
(and  rotational)  energy  is  then  conducted  to  the  tube  wall  with  a  decay  time,  x^c.  This  x^c  has 
been  calculated  according  to  the  equation 


ttc  = 


R  pCvK 
2.4\ 


(3.1) 

where  R  is  the  tube  radius,  Cy  is  specific  heat  at  constant  volume,  p  is  gas  density  and  A  is  the 
thermal  conductivity,  and  K  is  an  emperical  constant.  One  would  expect  xyc/xtc  to  equal  a 
constant,  but  until  now  little  has  been  determined  about  the  value  of  that  constant.  During  the 
past  year,  considerable  time  has  been  spent  in  trying  to  improve  the  method  of  fitting  this 
relaxation  model  to  the  experimental  curves  shown  in  Fig.  3.2  and  in  fitting  the  gains  calculated 
by  the  equations  of  reference  2  with  the  measured  gains.  These  efforts  include  the  following: 

1 .  The  temperature  dependence  of  the  vibrational  specific  heat  has  been  included  in 
the  relaxation  equations  (Eq.  2  of  Ref.  4). 

2.  The  two  channel  oscilloscope  has  been  replaced  with  a  four  channel  one  to  allow 
for  the  simultaneous  measurement  of  the  discharge  energy  and  acoustic  pulse. 
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3.  Changes  have  been  made  to  avoid  overdriving  the  microphone.  This  is  the  most 
important  change  to  date.  The  1  inch  B&K  microphone  that  was  in  use  last  year 
was  being  over  driven.  This  evidently  was  the  cause  for  some  of  the  very  large 
initial  gains  that  were  observed  following  the  electrical  discharge.  The  1  inch 
microphone  has  been  replaced  by  a  1/2  inch  microphone  and  care  has  been  taken 
to  make  sure  it  is  not  over  driven.  The  gains  being  measured  with  this 
microphone  are  still  larger  than  theory  but  are  now  of  the  same  order  of 
magnitude. 

4.  Previous  measurements  made  with  the  1  inch  microphone  have  been  repeated 
with  the  1/2  inch  microphone. 

5 .  Using  the  reasoning  involved  in  deriving  the  Euchen  expression, 
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it  has  been  shown  that 

TVc/xtc  =  1  +  2.25  R/C^  ,  (3.2) 

where  k  is  thermal  conductivity,  r|  is  viscosity,  cv  is  specific  heat  at  constant 

volume,  y  is  Cp/cv,  R  is  the  gas  constant  and  Coo  is  the  non-relaxing  part  of  the 

specific  heat  which  in  the  case  of  diatomic  gas  is  5/2  R.  This  is  an  important 
result  and  has  been  verified  by  the  curves  in  Fig.  3.2. 

The  effect  of  the  modulation  produced  by  the  sound  wave  upon  the  rate  at  which 
heat  is  conducted  to  the  tube  wall  has  been  included  in  the  equations  used  to  get 
the  sound  amplification.  This  does  not  significantly  effect  the  calculated  gains. 


Through  these  efforts  it  has  been  possible  to  fit  the  translational  temperature  curves  (Fig. 
3.2)  in  all  of  the  mixtures  with  a  single  adjustable  parameter  (the  K  in  Eq.  3.1  is  set  equal  to  0.7). 
Equation  3.1  is  derived  assuming  that  the  temperature  distribution  across  the  tube  varies  as  the  zero 
order  Bessel  function  with  the  first  root  at  the  tube  wall.  This  assumption  is  certainly  invalid  for 
times  soon  after  the  discharge.  The  0.7  value  for  K  is,  therefore,  reasonable.  The  fact  that  all  of 
the  data  can  be  fit  by  a  single  adjustable  parameter  is  especially  impressive  since  the  relaxation 
times  for  N2  in  the  various  mixtures  varies  by  a  factor  of  1400  and  the  rate  of  change  of  rv  with 
temperature  varies  by  two  orders  of  magnitude. 
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Having  established  the  value  of  K  in  Eq.  3.1  and  the  validity  of  Eq.  3.2,  it  is  now  possible 
to  extract  relaxation  times  as  a  function  of  temperature  from  the  curves  in  Fig.  3.2.  In  this  way  we 
have  determined  that  the  relaxation  time  for  N2  in  water  vapor  is  essentially  independent  of 
temperature  (a  point  previously  contested)  and  that  an  earlier  reported  room  temperature  value  is 
20%  too  large.  This  information  is  needed  in  predicting  sound  absorption  in  air  at  low  humidities 
and  high  temperatures. 

While  some  progress  has  been  made  in  getting  agreement  between  measured  and  calculated 
amplification  coefficients,  the  measured  values  are  still  larger  than  expected.  This  is  especially  true 
in  N2/H2O  mixtures.  Something  maybe  going  on  that  is  yet  not  understood.  Work  continues  in 
trying  to  understand  these  differences. 

3.4  Work  Planned  for  the  Coming  Year 

The  results  of  the  measurements  in  N2/He,  N2/CH4,  and  N2/H2O  should  be  submitted  for 
publication  before  the  end  of  the  year. 

In  order  to  predict  the  SACER  effect  it  is  necessary  to  know  the  vibrational  temperature.  At 
present,  as  explained  above,  this  vibrational  temperature  is  only  inferred  from  the  energy  in  the 
discharge  and  the  variation  of  the  translational  temperature.  It  will  greatly  aid  in  understanding  the 
experimental  results  if  this  vibrational  temperature  can  be  measured  directly.  Initial  efforts  have 
been  made  to  do  this  using  CO/H2  mixtures.  These  efforts  will  continue.  In  addition,  an  effort 
will  be  made  to  add  trace  amounts  of  an  infrared  active  gas  (perhaps  CO  or  CO2)  to  the  N2 
mixtures  that  are  strongly  UV  coupled  to  the  N2.  In  this  way  the  infrared  radiation  from  the  trace 
gas  will  follow  the  decay  of  the  N2  vibrational  temperature. 

While  these  experiments  were  begun  in  an  effort  to  observe  SACER,  the  phenomena  is 
now  well  enough  understood  that  it  can  be  used  to  measure  transition  rates  and  their  temperature 
dependencies  that  are  otherwise  hard  to  determine.  For  example,  there  is  considerable  uncertainty 
in  the  temperature  dependence  of  the  N2/H2O  relaxation  time.  This  information  is  needed  in  order 
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to  calculate  sound  absorption  in  both  N2  and  air  at  high  temperatures  and  low  humidities.  The 
measurements  recently  completed  should  provide  valuable  information  about  these  rates. 

As  time  permits,  an  effort  will  be  made  to  observe  SACER  in  polyatomic  gases.  In  this 
case,  the  sound  will  be  high  frequency  and  the  vibrational  excitation  will  be  produced  by  irradiating 
the  gas  with  a  high  intensity  pulse  of  laser  light.  This  experiment  would  extend  into  1989. 
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